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Generation of anomalously energetic suprathermal electrons was observed in simulation of a high- 
voltage dc discharge with electron emission from the cathode. An electron beam produced by the 
emission interacts with the nonuniform plasma in the discharge via a two-stream instability. Efficient 
energy transfer from the beam to the plasma electrons is ensured by the plasma nonuniformity. The 
electron beam excites plasma waves whose wavelength and phase speed gradually decrease towards 
anode. The short waves near the anode accelerate plasma bulk electrons to suprathermal energies. 

The sheath near the anode reflects some of the accelerated electrons back into the plasma. These 
electrons travel through the plasma, reflect near the cathode, and enter the accelerating area again 
but with a higher energy than before. Such particles are accelerated to energies much higher 
than after the first acceleration. This mechanism plays a role in explaining earlier experimental 
observations of energetic suprathermal electrons in similar discharges. 

PACS numbers: 52.35.Qz, 52.40.Mj, 52.65.-y, 52.77.-j 


The collisionless relaxation of high-current electron 
beams and plasma jets is important for many applica¬ 
tions: solar disruptions [i|, collective stopping! intense 
electron beams for inertial fusion applications^, Q, colli¬ 
sionless shocks and collision of flowing plasmas in astro¬ 
physics 0^: generation of suprathermal electrons 0]- 
Such beams are common in laboratory plasmas. Elec¬ 
trons emitted by electrodes surrounding or immersed in 
the plasma are accelerated by the sheath electric field 
and become the electron beams penetrating the plasma. 
In plasma applications where controlling the electron ve¬ 
locity distribution function (EVDF) is crucial 0 , these 
beams are an important factor capable of modifying the 
EVDE and affecting the discharge properties. 


Recently, Xu et al. [ll| and Chen and Funk [12| re¬ 
ported an EVDF measured in a dc-rf discharge with 
800 V dc voltage which has not only a peak at 800 
eV corresponding to the electrons emitted from the dc- 
biased electrode, but also a significant fraction of accel¬ 
erated electrons with energy up to several hundreds eV. 
In Refs. 11, [l^, the acceleration is explained invoking 


nonlinear wave-particle interactions 13-1^. It is sug¬ 
gested that the 800 eV beam excites fast long plasma 
waves which decay parametrically into ion acoustic waves 
and short plasma waves with much lower phase velocity, 
then the short waves accelerate slow electrons of the low- 
temperature plasma bulk. A similar mechanism may ex¬ 
plain enhanced enereetic tails in the energy spectrum of 
electrons in aurora 0, 2^ 211. Other mechanisms invoked 
to explain the EVDF structure include phase-bunching, 
kinematic effect of electrons being trapped between rf and 
dc potentials and released towards the rf electrode during 
short period of time when rf voltage is small compared 
to the dc voltage [22 1. 


In order to further understand the mechanism of ac¬ 


celeration in the experiment of Refs, [llj, ll2|, a beam- 


plasma system was simulated using a 1D3V particle-in¬ 
cell code EDIPIC 231. The EDIPIC was extensively 


benchmarked against available analytical predictions for 
two stream instability in linear and nonlinear regimes. 
Simulation results discussed below show that the accel¬ 
eration may be caused not only by the nonlinear effects, 
but by the effects related to the plasma nonuniformity 
as well. Note that profiling the background plasma den¬ 
sity is a method to control the beam energy deposition 
in plasma In particular, one can excite the in¬ 

stability and achieve beam energy deposition in desir¬ 
able regions while suppressing the instability in unde¬ 
sirable regmns by employing convective stabilization ef¬ 
fects H, 25|. In Refs. 0, 2^28| an experimental study 
and particle-in-cell simulations of interaction of a strong, 
warm electron beam with an inhomogeneous, bounded 
plasma showed that the trapped Langmuir waves can de¬ 
termine the frequencies and areas of localization of the 
wave field. In those studies the beam energy was rela¬ 
tively low (40 eV) and the acceleration of plasma elec¬ 
trons was not observed. 


The present paper discusses the mechanism of gener¬ 
ation of suprathermal electrons found in the simulation. 
In this process, the beam excites long plasma waves in the 
plasma body, the long plasma waves are converted into 
short plasma waves at the plasma periphery where the 
plasma density decreases, and the short waves accelerate 
bulk electrons. Some of these electrons return into the 
acceleration area again due to reflections from the sheath 
regions near the plasma boundaries. The second inter¬ 
action with the short plasma waves yields electrons with 
anomalously high energies. Repeating the acceleration 
for the same particle is more efficient than a one-time in¬ 
teraction between particles and plasma waves in density 
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FIG. 1: (a) Profiles of electron density (red) and electrostatic 
potential (blue) immediately before the electron emission at 
the cathode begins, t = 49 ns. (b) Amplitude of electric field 
oscillations vs time at point x = 12 mm (marked by arrow A 
in (a)). Arrows 1 and 2 in (b) mark times when snapshots 
shown in Figures [2] and O respectively, are obtained. 


gradient areas considered e.g. in 

The simulated collisionless argon plasma is bounded 
between anode (x=0) and cathode (x=40 mm). The 
cathode potential is zero, the anode potential is 800 V. 
The cathode can emit electrons with a constant flux 
and a temperature of 1 eV. The initial state of the 
system shown in Fig. [T] is the result of a preliminary 
simulation which starts with uniform plasma density 
no = 2 X 10^^ m^, the emission from the cathode turned 
off, the electron temperature = 2 eV, the ion tem¬ 
perature Ti^o = 0.03 eV. The preliminary simulation 
lasts for 4000 ns. During this time, wide areas with 
plasma density gradients form as shown in Fig. [T^. The 
time counter in the simulations with the beam starts at 
t = 0 and the electron emission starts at t = 50 ns. 
The emission current is 20.13 mA/cm^, the relative beam 
density a = Ub/no in the plasma density plateau area 
(11 mm < X < 27 mm) is a = 3.77 x 10““^. 

In the density plateau area, the plasma waves exited by 
the beam have the wavelength close to 2TrvbjuJo^pe where 
<^o,pe = e{nojmeSo)^^'^ and the amplitude growing along 
the beam propagation towards anode [1^, see Figs, 
and Hi. In the density gradient area, a; < 12 mm, the 
profile of the electric field drastically changes: the ampli¬ 
tude and the wavelength decrease towards anode. Phase 
plots in Figs.H andH show that bulk plasma and beam 
electrons perform strong oscillations in the electric field 
of the plasma waves. Suprathermal electrons (STEs) are 
clearly seen in the EVDF near the anode, see Fig. |4l 
The electrons are accelerated from the initial Maxwellian 
EVDF with the electron temperature Tg.o = 2 eV to en¬ 
ergies up to 60 eV. 

The process of acceleration is demonstrated in Figs. [S] 
and M where the color map of the electric field E(t,x) 
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FIG. 2: (a) Profile of the electric field, (b) Phase plane 
“coordinate-velocity” of bulk electrons (colormap, left ver¬ 
tical coordinate axis) and the electron beam (blue markers, 
right vertical coordinate axis). The snapshots are obtained at 
time 114.97 ns marked by arrow 1 in Figure [T]3. 


is presented along with the trajectories of few electrons. 
The change of the slope of field pattern seen in the color 
map (Figs. H and H) corresponds to the decrease of 
the phase velocity towards the anode due to the plasma 
density inhomogeneity. The frequency of the wave ex¬ 
cited by the beam is the same everywhere, including the 
density gradient regions, and is very close to the plasma 
frequency in the density plateau oJo^pe- Assuming that 
the plasma wave satisfies the dispersion relation = 
Wpg -I- ik’^Te/m.e and setting k = uiQ^pe/vb in the density 
plateau, the wave frequency is a;o,pe(l + 3Te/TOe'i'b)^^^- 
Substituting this frequency into the dispersion equation 
and using /Uq, one can obtain the phase 

velocity of the wave as a function of coordinate 

s/iTe/rUe 


Vphix) = 


^1 - (1 - 3Telmevl)n{x)fno ’ 


( 1 ) 


where Tg and me are the electron temperature and mass, 
Ueix) is the plasma density which depends on coordinate 
X, and Vb is the beam velocity. The analytical expres¬ 
sion for the phase velocity o is in good agreement with 
the results of the numerical simulation, compare the blue 
curve with the red markers in Fig. H- 

After the emission from the cathode started, STEs 
moving towards the anode appear first downstream of 
the area with the strongest oscillations, see Fig. H for 
cc < 8 mm. At the peak of the instability (arrow 1 in 
Fig.dlD), maximal energy of these electrons reaches about 
34 eV, see curve 1 in Fig. SI While the accelerated elec¬ 
trons escape to the anode, the plasma potential relative 
to the anode increases. Some STEs are reflected by the 
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FIG. 3: Same as in Figure [2l but the snapshots are obtained 
at time 149.96 ns marked by arrow 2 in Figure [TJs. 



FIG. 4: Bulk electron velocity distribution function obtained 
using particles with 1 mm < x < 7 mm, the horizontal axis 
is in energy units, negative values correspond to propagation 
in the negative ^-direction. The red and blue curves are for 
the system state shown in Figures [2] and [3l respectively. 


FIG. 5: Interaction of test particles of different initial energy 
with the waves, (a) Color map of the electric field as a func¬ 
tion of coordinate and time Ex{x,t). The positive electric 
field shown with the red color is directed rightward, towards 
the cathode. The curves are the test particles trajectories 
x{t). (b) Energy of test particles vs coordinate. In (a) and 
(b), the black solid curves correspond to a particle with ini¬ 
tial energy 7 eV at a; = 20 mm at t = 110.77 ns, the cyan 
curves correspond to a particle with initial energy 4 eV at 
a; = 20 mm at t = 108.08 ns. The solid blue curve in (b) 
marks the electron energy corresponding to the theoretical 
local phase velocity of the wave given by Eq. [T] In (b), the 
negative values of energy correspond to propagation in the 
negative x direction (towards the anode), the arrow points in 
the direction of propagation of the test particles. 


enhanced potential barrier of the anode sheath and, after 
a round trip through the plasma with another reflection 
near the cathode, they return into the area with intense 
plasma oscillations. When it happens, the energies of 
STEs moving towards the anode downstream of the in¬ 
stability maximum increase even more, see Fig. [5 )d and 
curve 2 in Fig. |4l The maximal energy of STEs (reg¬ 
istered at t = 139.97 ns) is about 80 eV. At the first 
minimum of the electric field amplitude as a function of 
time (arrow 2 in Fig. [Ib), the maximal energy is lower, 
about 65 eV, but still almost twice the maximal energy 
of the one-stage acceleration, compare curves 2 and 1 in 
Fig.m 

Acceleration of bulk electrons by the short waves is 
studied with test particles which move in the electric 
field provided by the simulation. A test particle with 
low initial energy (4 eV) does not get net acceleration, 
its energy only oscillates along the particle trajectory, 
see cyan curves in Fig. [5l A particle with higher initial 
energy (7 eV), however, is accelerated by these waves, see 
the black curves in Fig.[S] Note that before and after the 


acceleration, the energy of the 7 eV test particle in the 
wave field oscillates, see the black curve in Fig. [SJd. The 
acceleration occurs if at the moment when the energy 
of the particle is maximal, the velocity of the particle is 
equal to the local phase speed of the wave, compare the 
black energy-vs-coordinate curve of the 7 eV test parti¬ 
cle and the blue curve corresponding to the wave phase 
velocity in Fig. [5 }d. 

The average energy of the 7 eV test particle after the 
acceleration is about 34 eV which is sufficient to pene¬ 
trate through the sheath potential barrier (about 23eV). 
Another test particle with initial energy of 6 eV after the 
first acceleration (point A in Figs. and [B}:) has only 
20 eV. This particle is reflected by the anode sheath, 
travels through the whole plasma, reflects near the cath¬ 
ode, and approaches the instability area with the energy 
of about 20 eV, see Fig. |Bb. The particle is accelerated 
again at the segment between 9.5 mm and 11.5 mm, see 
point B in Figs. Eb and Et. The energy of the particle 
after the acceleration is about 65 eV, which is more than 
ten times the initial particle energy. The first accelera¬ 
tion occurs similar to the case described in the previous 



















4 


E^(10'‘v/m) 


-10 -8 -6 -4 -2 0 2 4 6 8 10 



6 7 8 9 10 11 12 13 14 



X (mm) 


FIG. 6: Two-stage acceleration of a test particle. The initial 
energy is 6 eV at x = 20 mm at t = 110.44 ns. (a) The 
first stage: color map of the electric field as a function of 
coordinate and time E^ix, t). The positive electric field shown 
with the red color is directed rightward, towards cathode. 
The curve is the test particle trajectory x{t). (b) Same as (a) 
but for the second stage, (c) Energy of the test particle vs 
coordinate. In (b), the negative values of energy correspond to 
propagation in the negative x direction (towards the anode). 
Point O marks the starting point, point A - acceleration in 
the first stage, point B - acceleration in the second stage, and 
point E - escape of particle at the anode. 


paragraph except that due to the lower initial particle 
speed the acceleration occurs closer to the anode where 
the wave amplitude is lower, compare Figs.lBjD and[S^. 

The second acceleration occurs in the area immedi¬ 
ately adjacent to the density plateau where the wave 
phase speed and, in general, amplitude are higher. It 
is necessary to mention, however, that by this time the 
short waves in the majority of the density slope area de¬ 
cayed because they were involved in intense acceleration 
of plasma particles, see Figs. [5)3 and [5^ for a; < 9 mm, 
and the instability amplitude is also noticeably lower 
than during the first acceleration, compare wave fields 
in Figs. and HJa. The decay of the plasma waves is 
a significant limiting factor for the observed two-stage 
accelerating mechanism. 

Summarizing, in a dc plasma-beam system with 
nonuniform density, suprathermal electrons are gener¬ 
ated by short plasma waves excited at the density gra¬ 
dients. Some of the accelerated electrons may be re¬ 


flected by the anode sheath and reintroduced into the 
two-stream instability area where they will be acceler¬ 
ated one more time. The energy of an electron after the 
second acceleration can be an order of magnitude higher 
than its initial energy and few times more than its energy 
after the first acceleration sta ge. The simulation above, 
unlike the experiments of ll|, ll2| , does not include the 
rf voltage. This has both favorable and negative conse¬ 
quences. On one hand, variation of the energy of the 
beam by the rf voltage may reduce the efficiency of exci¬ 
tation of plasma waves. On the other hand, in a real dc-rf 
system, plasma bulk electrons are trapped (at least part 
of the rf period) in a potential well much deeper than the 
one created by thermal electron motion. This will pre¬ 
vent suprathermal electrons from escaping to the walls 
and allow them to go through two or more acceleration 
cycles as described above. 
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